Preparing and probing atomic number states with an atom interferometer 
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We describe the controlled loading and measurement of number-squeezed states and Poisson 
states of atoms in individual sites of a double well optical lattice. These states are input to an atom 
interferometer that is realized by symmetrically splitting individual lattice sites into double wells, 
allowing atoms in individual sites to evolve independently. The two paths then interfere, creating a 
matter-wave double-slit diffraction pattern. The time evolution of the double-slit diffraction pattern 
is used to measure the number statistics of the input state. The flexibility of our double well lattice 
provides a means to detect the presence of empty lattice sites, an important and so far unmeasured 
factor in determining the purity of a Mott state. 

PACS numbers: 03.75.Gg, 03.67.-a, 32.80.Pj 



The optical beam splitter, with its two input and two 
output modes, is one of the simplest examples of a two 
mode quantum system. At the quantum level this funda- 
mental system has interesting, nonclassical behavior such 
as the quantum interference between correlated, indistin- 
guishable photons This system becomes even richer 
when the particles interact. While photons may effec- 
tively interact in nonlinear media, the atom optics ana- 
log is naturally interacting. The two mode beam split- 
ter has already been applied in atom optics experiments; 
several experiments have split a trapped Bose-Einstein 
condensate (BEC) by raising a barrier to separate the 
condensate into two independent condensates [1, H, 0> S] • 
In these experiments the number of atoms is large, in 
a regime where few-particle quantum interference effects 
cannot be seen. Here we demonstrate a new few-atom 
quantum "beam splitter" [f| and use it to create and ana- 
lyze classical and nonclassical states, resulting in interest- 
ing few-atom quantum effects. The ability to create and 
analyze such states provides a probe of many-body states 
in a lattice, a platform for fundamental studies of few- 
particle, interacting systems, and is of paramount impor- 
tance for quantum computation with neutral atoms. 

We realize an atomic analog of the optical two-mode 
quantum beam splitter with 87 Rb atoms loaded into an 
optical lattice of double wells [tJ- This 3D lattice has a 
unit cell that can be dynamically transformed between 
single well and double well configurations (Fig. QJ,). In 
analogy with the optical beam splitter, the input and 
output modes in this lattice are either the ground \g) 
and first excited |e) state of the single well or the ground 
state of the left \L) and right \R) sites of the double well. 
Dynamically switching between these two configurations 
effectively creates an array of "beam splitters" (BS) cou- 
pling these modes (in the limit that higher modes are not 
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FIG. 1: (color online) (a) Action of a double well beam split- 
ter. Different colors indicate the lack of phase coherence be- 
tween double well sites, (b) Optical analog of our atom inter- 
ferometer with an example diffraction pattern resulting from 
a single input mode. 



excited). The interactions between atoms and our con- 
trol of the beam splitter time scale extends the possible 
actions beyond those of a traditional beam splitter. 

The speed at which the beam splitter is applied to the 
atoms determines whether interactions play a role in the 
splitting: if the BS is applied quickly compared to inter- 
action energies, the interactions are unimportant. This 
"non-interacting" beam splitter (NI-BS) is analogous to 
the linear optics case: atoms in a single input mode (for 
example, |<7)) are divided into two output modes (here, 
\L) and \R)) according to binomial statistics. In con- 
trast if the BS is applied slowly compared to interac- 
tion energies, the interactions change the statistics of the 
distribution between output modes. This "interacting" 
beam splitter (I-BS) has no simple optical analog. We 
use the NI-BS and I-BS to create and analyze classi- 
cal Poisson states and interesting nonclassical few-atom 
squeezed states. The ability to controllably and reliably 
load high fidelity Fock states is attractive for applications 
in quantum information, including the implementation of 
a two-qubit operation in an optical lattice. 
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FIG. 2: (color online) (a) Dumping the population initially 
in \L) into highly excited, possibly unbound bands, (b) The 
resulting TOF image: the population initially in \R) appears 
in the first BZ of the A- lattice; the population initially in \L) 
appears in higher BZs as a ring around the first BZ. (c) Cali- 
bration of the BS in Fig. [T] After applying the BS at a given 
value of Sx/X, the population in the \L) (red) or \R) (blue) 
sites are measured by dumping the other site and measuring 
the fraction of atoms remaining in the first BZ. The value of 
Sx/X where the two curves cross corresponds to a 50/50 BS. 



We probe these states using atom interferometry 
(Fig. |T|) . The NI-BS symmetrically and coherently splits 
a single lattice site into a double well. Atoms in indi- 
vidual sites of a double well evolve independently, anal- 
ogous to the individual arms of an interferometer. We 
interfere the two paths by releasing the atoms and al- 
lowing them to expand and overlap during time-of-flight 
(TOF) The interference from a single pair of sites results 
in a double-slit diffraction pattern. In this experiment 
the initial states were prepared to have no phase coher- 
ence between unit cells; thus the observed image is an 
incoherent sum of doublc-slit interference patterns from 
each of the double wells. 

The input states were created and analyzed in a dou- 
ble well optical lattice Q. We applied this lattice to a 
BEC of 87 Rb atoms in 55 1/2 \F = l,m F = -1} created 
as in Q. The lattice comprises two independent 2D lat- 
tices having periodicity along x of either A or A/2 (Fig. 
[TJl) where A = 815 nm which is red-detuned from the 
87 Rb D-lines. The relative position Sx of the "A-lattice" 
and the "A/2-lattice" and their intensities can be dy- 
namically controlled. Transforming from the A-lattice 
to a predominantly A/2-lattice with the minima of the 
A/2-lattice symmetrically straddling the minima of the 
A-lattice accomplishes the splitting shown in Fig. Q] We 
use an independent ID "vertical lattice" in the z direc- 
tion to confine atoms to sites of a 3D lattice. The double 
well and vertical lattices are focused to 1/e 2 Gaussian 
beam radii of 170 /im and 250 /im respectively and have 
a relative detuning of 160 MHz. The final depths of the 
A/2-lattice and vertical lattice are w 30Er and w 40£'fl, 
respectively (En = h 2 k 2 /(2m) — h x 3.5 kHz, k — 2tt/\, 
m is the Rb mass) . The atomic momentum distribution 
in the x, y-plane is absorption-imaged after TOF. 

We measure the populations in the two output ports, 
\L) and \R), by selectively imparting energy to atoms in 
one of the two sites, which separates the populations in 



TOF (Fig. HJ). Transforming to the A-lattice (in 300 (is) 
with its minimum centered over the \R) site, imparts en- 
ergy to \L) atoms; \R) atoms remain in the ground state 
of the A-lattice (Fig. [5^). The final states' Brillouin zones 
(BZs), which are filled, are observed (see for instance 
0j@]) by turning off the lattice in 500 /is and imaging af- 
ter TOF (Fig. [2b>). The atoms in the center of the image 
correspond to the filled first BZ of the A-lattice 0] and 
were originally in \R). The atoms in the ring correspond 
to higher BZs and were originally in \L). This method 
allows for direct measurement of the atom population in 
\L) and \R). Using this method we can accurately cali- 
brate the relative position of the two lattices that gives 
equal splitting of the NI-BS (shown in Fig. [3b). This 
50/50 BS was used for all experiments described here. 

To prepare number squeezed states with N = 1 atom 
in \g) of the A-lattice sites, we used a two-stage "slow" 
loading procedure. We started with w 2 x 10 4 atoms 
from a BEC with no discernable uncondensed fraction in 
a magnetic trap with loj_/2it = 24 Hz and u)\\/2tt = 8 Hz 
and Thomas Fermi (TF) radii of w 13 /im and ks 32 /im, 
respectively. The vertical lattice was slowly raised in 
250 ms, forming an array of 2D "pancakes." After the 
first 100 ms of this loading, we raised the A-lattice in the 
remaining 150 ms [la ]. Ideally, at zero temperature, suffi- 
ciently slow loading produces a vertical array of 2D Mott 
insulating systems with one atom in \g) of each A-site. To 
prepare Poisson states with (N) « 1 in \g), we used the 
same procedure, except we loaded both lattices quickly 
(500 /is). This "fast" load is adiabatic only compared to 
band excitation [7|. We destroyed the phase coherence 
between atoms in different A-sites by holding atoms in 
the lattice for 5 ms. After preparing either of these ini- 
tial states, we applied the NI-BS by transforming to the 
A/2-lattice in 300 /is, a timescale fast enough to maintain 
coherence within a double well, yet slow enough to avoid 
vibrational excitation. After a variable hold time t, the 
atoms were released and imaged after 13 ms TOF. 

Fig. [3bo shows an example doublc-slit interference pat- 
tern. We extract the fringe visibility C(t) and the spatial 
phase 4>(t) by fitting the profile to: 



F(x,t) = Ae 2, 
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where a is determined by the width of the ground state 
wavefunction in a A-site, and 27rA is the fringe spacing. 
We use C (t) to measure the number statistics of our input 
states. 4>(t) evolves linearly due to an e nerg y offset (tilt) 
V between the sites within a double well [l6| (see Fig. [3)}). 

The Nl-beam splitter binomially splits N atoms in each 
A-lattice site into superpositions of states \ni,,nR), with 
til atoms in \L) and ur = N — til atoms in \R), resulting 
in the time dependent wavefunction 
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FIG. 3: (color online) (a) A double-slit interference pattern 
(inset) and integrated density profile after 13 ms TOF. (b) 
Example y-integrated diffraction pattern as a function of time 
for tilts V = 4.90(3) kHz (left) and V = -4.90(3) kHz (right), 
(c) C(t) for N = 1 (black o), N = 2 (green x), and (N) = 1 
(blue ■). Solid lines are a fit to Eq. H Note: The plotted 
C{i) for N = 2 was scaled by 2. (d) 0(t) for the cases shown 
in (c). The <f>(t) points which were derived from a signal with 
2 X C(t) < 0.1 are shaded in light green. <f>(t = 0) = for all 
three cases, but each has been shifted for visual clarity. 
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where \c(nL,n R )\ — (h) {N\/riL\n R \). In the two- 
mode approximation interactions and tilt determine 
hw(n Ll n R ) = n L V+{l/2)Un L (n L - 1) + (l/2)Un R (n R - 
1), where U is the on-site interaction energy per particle. 
For our system U/h w 3 kHz. 

In the absence of interactions and inhomogeneities, 
\tp) always factors into products of single particle states, 



IV> = (|i,o> + . 



0, 1))^, where uj tm = V/h. All the 



terms in Eq. [T] constructively interfere and C (t) is max- 
imized and constant. When interactions are included, 
is not always factorable, and C(t) (for atoms in the 
ground band) is modulated by ^os 1 * -1 (Ut/H)\ [13, 11 1. 
Summed over a distribution of site occupation numbers, 
the total C{t) is given by 



-rt 



C(t) = 



N„ 



(ut/h) 



N=l 



(2) 



where Jn is the fraction of A-sites with N atoms (for 
our case we limited N mgx — 4), /jv = 1, and T is 
an empirical dephasing rate that we attribute to lattice 
inhomogeneities. Fitting C(t) with Eq. [2 we can extract 
the relative fractions /jv- 

For the slow loading case with Nbec < 2 x 10 4 (aver- 
age filling < 1 per site) we expect to have a Mott state; 
we measure C(t) to be approximately constant (see Fig. 
|3J;), indicating that occupied A-lattice sites have N = 1. 



(C(t) shows a slow decay, which we attribute to tilt in- 
homogeneities.) Extracting /jy from C(t), we determine 
fx = 0.94(6). (As shown in Fig.0^,, fi was optimized at a 
load time of 150 ms.) The uncertainty in fi is dominated 
by the shot-to-shot scatter in the double-slit visibility C, 
which is significantly larger than for fast loading. We do 
not understand this increased noise, but note that the 
shot-to-shot scatter in <f) (Fig.EJi) is not increased. This 
suggests an increased sensitivity to an uncontrolled ini- 
tial condition, e.g., number or temperature. 

The near-zero value of /2 indicates a strongly number- 
squeezed initial state, as expected in a Mott insulator 
(compared to the Poisson value of 0.3 for average fill- 
ing of unity, ignoring f ). The value of /1 does not, 
however, necessarily represent the fidelity of the Mott 
state since this measurement is insensitive to sites with 
N = 0. Assuming fx = 0.94 and fa = 0.06 and a ho- 
mogeneous system, we can roughly estimate the temper- 
ature T in the A-lattice from fa = e -U/k B T _ gives 
T < 0MU/k B ^ 50 nK. 

For the fast loading case where we expect a Poisson 
distribution of site occupation numbers, there can be 
N > 1 atoms at any A-site. In this case we find that 
C(t) shows distinct collapses and revivals, caused by in- 
teractions (see Fig. |3j:). Only at integer multiples I of 
T rev = h/U does the wavefunction within a site factor: 

1]) = ( 1 1 , 0) + e Uu)tmTx " v |0, 1}) . We observe revivals in 
C(t) with a period of ~ 350 [is, in good agreement with 
our calculated value oiU/h = 2.88 kHz. While similar to 
the interaction-induced collapse and revival seen in [13], 
the interference here is between the sites of a double well 
and is created by a topology change in the lattice, rather 
than between sites extending over the entire lattice and 
created by changing the lattice depth. This allows us to 
see interference even in the absence of global coherence, 
as evidenced by the interference pattern generated from 
the highly squeezed N = 1 slow load state. 

Using Eq. [5] to fit the fast load C(t), we extract the 
number distribution of atoms in A-lattice sites. By fit- 
ting that distribution to an average over Poisson distribu- 
tions whose local mean occupation numbers are propor- 
tional to the initial TF-density profile, we extract (N), 
the peak number density of atoms in A-lattice sites. This 
weighted Poissonian distribution fits the data within ex- 
perimental uncertainty. Fig. [4]b shows (N) as a function 
of A^bec and the power law scaling (N^Q C ) for (N) ex- 
pected from the TF-approximation. For large values of 
A<bec the data deviates from the trend set by the lower 
values of A^bec- While not completely understood, this 
result may be attributed to three-body loss mechanisms 
and a larger thermal fraction for larger values of A^bec- 

We now consider input states with N = 2 atoms in \g) 
of the A-lattice sites. If we simply loaded the A-latticc 
from our harmonic trap with increased initial peak den- 
sity, the resulting "shell structure" should produce a sig- 
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the wavefunction factored into oppositely-phased single 
particle states, \tp) = (|1,0) - e ^t ilt T rov /2 |n 1)) 2 ; giving 
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FIG. 4: We use Eq. [2] to fit C(t) for the fast and slow load 
data, (a) Extracted /i for varied slow load times. The uncer- 
tainties are from the fit to Eq. [2] (b) Using Eq. [2] to fit the 
fast load C (t) we extract values for /i , /2 , fs, and f$ . From 
those coefficients we determine {N ) , the peak number density 
of atoms in A-lattice sites, as a function of Abec- The line 
and the grey shaded area are the peak number density scal- 
ings predicted by the TF-approximation (with no adjustable 
parameters) and our estimated uncertainty in the calibration 
of Nbec (factor of ±1.5). The uncertainties on the individual 
points are from a combination of the uncertainties in fi and 
from the fit determining (N). 



nificant number of sites with N — 1 surrounding the cen- 
tral core of N — 2 [3, [3] , limiting the maximum value 
of f-2- In principle, the flexibility of the double well lattice 
can be used to construct a purer N — 2 state in the A- 
lattice by combining two A/2-lattice sites: we first load 
the vertical lattice as in the slow load case above, but 
with an initial number (ATbec ~ 7 x 10 4 ) large enough 
to produce ~ 1 atom per site in the A/2-lattice. We 
then raise the A/2-lattice in 100 ms (ideally, slow enough 
to make a Mott insulating state with one atom in the 
ground state of each A/2-site). Finally we adiabatically 
(in 30 ms) combine pairs of A/2-sites into A-sites, thus 
performing an I-BS. The input states are twin Fock states 
of one atom each in \L) and \R), the two-atom interacting 
ground state of the A/2-lattice. The output modes are 
\g) and |e) states of A-sites. Interactions and adiabaticity 
ensure that both atoms go into \g), the two-atom ground 
state of the A-lattice. This "constructed pair" state is 
then input to our interferometer (split with an NI-BS). 

As shown in Fig. [3]:,d we find that after using the con- 
structed pair technique to load the A-lattice, C{t) and 
4>{t) show unique features which can only be attributed 
to A-sites with N — 2. As in the fast load case, C(t) 
collapses and revives, but there are revivals at half the 
period compared to the fast load case (T rcv /2 rather than 
T rcv ). Like the slow and fast load cases, <f>(t) evolves 
linearly as wtiit (here ps 0), but here <j>(t) makes clear 
jumps between revivals. These jumps are ~ n, partic- 
ularly if only the <p(t) points with 2 x C(t) > 0.1 are 
considered. At lT rcv the wavefunction can be factored as 



a 7T phase-shifted interference pattern 13, U | ■ 

Our maximum measured C in Fig. [3b for the con- 
structed pairs technique is a factor of ~ 2 less than the 
fast and slow load cases, and our measured C(T rev /2) 
is lower than C(T rcv ). These differences can be largely 
explained by the fact that after the I-BS, a large frac- 
tion of atoms (~ 30%) are in |e) (revealed by imaging 
the A-lattice BZ), and make an interference pattern n- 
out-of-phase with that made by atoms in \g) [17|. After 
the NI-BS the incoherent sum of sites with atoms equally 
in \g) and |e) — in total ~ 60% of the atoms — thus con- 
tributes C(t) = 0. |e)-population could be caused by 
I-BS input states such as \riL,nB,) = |2,0), |2,1), and 
1 1 , 0) , all of which would reflect impurities in the initial 
A/2-lattice Mott state. The presence of atoms in |e), 
however, also points to a unique capability of the double 
well lattice: traditionally it has been impossible to detect 
empty sites, but by using the constructed pairs technique 
to load the A-lattice and measuring N = 1, we can detect 
the presence of empty sites in the initial A/2-lattice. 

Using the versatility and control of the double well lat- 
tice, future studies can also be focused on creating and 
exploring more complicated input states with interferom- 
etry, such as entangled |iV, 0) + |0,7V) states. Using the 
constructed pairs technique but with a NI-BS to com- 
bine pairs of A/2-sites to single A-sites, we can create 
|2,0) + |0, 2) in the \g,e) basis. This would be the atom 
optics equivalent of the Hong-Ou-Mandel effect [l[ . 
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|0, 1)) , so C(t) revives. Unlike 



the fast load case, at ^T rov /2, C(t) also revives but with 
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[15] We load to ~ 12 E R in 50 ms then to ~ 100 E R in 100 ms. 

[16] V can result from misalignment of the A/2-lattice 

which we can in principle compensate with a small 
amount of the A-lattice. The residual ~ 1 kHz tilt in 
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Fig. [3}! had no effect on the results presented here. 
[17] Eq. [5]does not account for [e) population. 



